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ABSTRACT 


Conventional near-shore large-scale hydrographic surveys 
use only two lines of position (LOPs) for position fixing. 
Previous works have proved that accuracies needed to meet 
the present International Hydrographic Organization 
standards are frequently not achieved using conventional 
surveying methods. The concept of using multiple lines of 
position (MLOP) adjusted by the least squares method was 
described. Actual field measurements acquired in the autumn 
of 1984 were processed to ascertain the increase in accuracy 
using MLOP versus conventional two LOPs on each hydrographic 
position. Recommendations a to use a four-range fully 
automated position-fixing method to increase production, 
improve data quality, and have better control of the survey 
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I. INTRODUCTION 


A. BACKGROUND 

À hydrographic survey is a compilation of numerous types 
of data which are acquired to ultimately produce a nautical 
chart. The two most important нр. of a hydrographic 
Survey are (1) the depth sounding of the seabed and (2) the 
Simultaneous acquisition of a horizontal position which is 
attached to the sounding. If either of these two aspects 
are in error then the overall accuracy of the survey is 
Posts This thesis will address the latter aspect.  Ingham 
[1984, p. 14] states 

A major element of the hydrographic survey is the . 
dynamic positional control of the survey vessel. 

Indeed, the great majority of the surveyor's tasks 
offshore today require positioning services alone. 

The accurate determination of a position at sea has 
always been a painstaking effort for the hydrographic 
Surveyor. First, a survey of the coast adjoining the 
offshore survey area must be made to establish an accurate 
position of each shore control site. This shore control 
must then be "extended" offshore by a variety of means and 
systems depending on the scale and purpose of the survey. 
The position of a survey vessel is either determined Ый а 


visual method, such as a three-point sextant fix to visual 
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Signals onshore, or by a radio navigation method. Іп the 
very near future, the Global Positioning System (GPS) 
promises to completely revise the method of offshore 
navigation and surveying. Consistent differential 
positional accuracies of about 10 m, or less, are expected 
in real time once the system becomes fully operational thus 
allowing the hydrographer the freedom of surveying as far 
offshore as the situation demands. 

This thesis will address various ways to increase the 
positional accuracy of nearshore Surveys. Nearshore areas 
(i.e., channels, harbors, and bays) pose the most hazard to 
the mariner and are the areas where accurate position 
determination of the bottom features is essential. 

ЕЕ ана лбу of 1:5,000 scale, and larger, require 
intensive and time consuming delineation of shoreline 
features and bottom topography with positional accuracies of 
5 т, or less. Presently, GPS will not provide these 
accuracies in a real-time mode. Conventional hydrographic 
Survey methods are the only way to determine the position of 
a survey vessel in real time. Mobley stated [Saxena, 1980, 
p. A-4], 

Positioning systems for precise navigation for survey 

operations are barely adequate for the scales of surveys 

being accomplished today. A lot of time and effort is 

lost due to equipment down time, moving equipment and 

the use of only 2 lines of position . . . which the 

present National Ocean Survey (NOS) DAS can only accept. 


New and future replacement systems providing redundant 
LOPs are a must if higher accuracies are desired 
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offshore. Global Navigation may be the answer but will 
be 5 to 10 years before readily available to all users. 

If conventional systems and methods are to be used, a 
redundancy in the measurements 15 necessary to obtain а 
higher order of accuracy than the present system of 
acquiring two lines of position (LOPs) on each fix. This 
thesis will investigate the application of multiple lines of 
position (MLOP), adjusted by the least squares method, for 


hydrographic position control. 


B. HYDROGRAPHIC STANDARDS OF ACCURACY AND METHODS 
The National Ocean Service (NOS) is the agency 
responsible for nautical charting in the United States and 
has documented many standards and procedures for the 
acquisition of data during a hydrographic survey. With 
respect to the position of the survey vessel, the 
Hydrographic Manual [Umbach, 1976, p. 4-14] states 
Hydrographic position fixes on sounding lines are almost 
always determined by the intersection of two lines of 
position; dead-reckoning positions based on course and 
speed provide an additional internal check. 
Also Ingham [1984, p. 15] states in his chapter on 
Position-fixing Afloat 
For any fix one requires at least two lines of position 
AE obtained from readily identifiable sources. These 
should intersect at the largest angle possible for good 
fix-geometry, and have the smallest possible standard 


deviation, resulting in the smallest area of 
uncertainty. 
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NOS requires that for surveys of 1:5,000 scale, and 
larger, at least one of the LOPs must be from a visual 
method. Thus, a standard way of acquiring positional fixes 
on large-scale surveys has been the range-azimuth 
(rho-theta) method as addressed by Waltz [1983]. One of the 
reasons this hybrid method is popular is because the two 
LOPs will always intersect at 90°, the best intersection 
possible for good accuracy. 

Unlike mariners, the hydrographer has access to more 
sophisticated and more accurate positioning systems. But as 
Riemersma [1980, p. 10] observes 

There is no mariner or surveyor in the world who in 

conventional navigation does not use a check bearing 

just to make sure. Why then does one, with 

radiopositioning [sic] systems, use only two LOPs? 
By using only two LOPs there exists a level of uncertainty 
with each position which can only be resolved during the 
post processing phase of the survey [Perugini, 1984]. In 
the past, two LOPs were the most feasible and sometimes the 
only way to determine the position of the survey vessel. 
With the rapid advance in electronic technology, the 
acquisition of multiple-range LOPs are as common today as 
two-range LOPs were 20 years ago. However, due to cost 
limitations, or possibly the aversion of changing tried and 
true operational procedures, hydrographic surveys are still 


usually controlled with only two LOPs. The acquisition of 
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only two LOPs on each fix will change in the very near 
future for NOS as they are in the process of implementing a 
specification from the 1982 International Hydrographic 
Organization (THO) (Anonymous, 1932) p. 5] stating 

It is desireable that whenever positions are determined 

by the intersection of lines of position, three such 
` lines be used. The angle between any pair shoaid noe. 

Less thangs 0 ae 

With MLOP the hydrographer will have redundant 

observations on each fix. Босс redundant measurements may 
be adjusted, as in normal land surveying, by the least 
squares method to acquire the "best" fit. Applying the 
least squares method to hydrographic position adjustment 
allows the hydrographer to statistically determine errors in 
real time if computer programs are operated aboard the 


survey vessel [Silva, 1982, p. 100]. 


С. “ОБЛ 5 

Will the use of MLOP increase the accuracy of the fix 
over the present method of two LOPs? This question has been 
addressed by both Kaplan [1980] and Silva [1982] who found 
that MLOP do improve the accuracy of the survey. However, 
the method was applied to theoretical data only. Therefore, 
one objective of this thesis was to acquire data for a 
portion of an actual hydrographic survey with a minimum of 


four LOPs on each position fix. A model was constructed to 
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determine the least squares position for each fix from the 
MLOP. The precision of each MLOP fix was compared to the 

precision of a fix obtained by two LOPs to determine if a 

major difference did exist between the two methods. 

To obtain the additional LOPs, a significant increase in 
the cost of operations can be expected and must be 
justified. More LOPs will require more horizontal control 
for shore reference points. Additional survey equipment, 
electronic range measurement units or azimuth measurement 
units, or both, will be necessary to acquire the four LOPs. 
Initial equipment costs could possibly double at the onset 
of this method of surveying. Consequently, all of the costs 
associated with the MLOP method must be weighted against the 
optimal position accuracy requirement to determine if the 
additional accuracy is feasible within the scope of a normal 
hydrographic survey. ` 

A final objective of this thesis was to determine an 
optimal system which is the most cost effective as well as 
one which will meet the new IHO standards for position 


fixing of soundings. 
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II. PROBLEM DEFINITION 


The geometry of a standard position fix using two LOPs 
must be understood before the various MLOP methods may be 
evaluated. Therefore, the configuration of large-scale 
positioning methods, the errors unique to the equipment 
used, and the propagation of these errors into the offshore 


position warrant discussion in this chapter. 


A. LARGE-SCALE SURVEY POSITIONING METHODS 

The intersection of two lines defines a point in space.- 
In navigation and offshore surveying each of the two lines, 
or LOPs, can be generated via an electronic or visual 
method. Traditional hydrography requires that the two LOPs 
generated intersect at an angle greater than 30? and less 
һап 150°. | 

ре; used for computations are: 


a) X coordinate positive to the east, Y coordinate 
positive to the north 


b) Occupied Shore Stations; numerals l, 2, 


c) Unoccupied Shore Station used for Initial Azimuths; 
numeral Q 


d) Inverse Distance between shore stations; D 
e) Ranges to vessel, point P; Ry, Ro, 

Е) Computed horizontal angle; aq 

g) Measured horizontal angles; a); ag; 

h) Azimuths of lines; 01.2» О1,Р» 
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l. Range-Range Positioning Method 


An electronic range-range positioning method on 
large-scale surveys is generally accomplished with 
short-range (line-of-sight) equipment. At these short 
distances the ranges to the vessel are determined by 
measuring the time of transit or difference in transit times 
of an electromagnetic pulse [Davis et al., 1981]. The 
transit-time measurement starts with an interrogation pulse 
from.a master transponder onboard the vessel which is 
received by a remote transponder erected over a known 
control point and returned to the master transponder. The 
total transit time is then converted to a distance by 
assuming an average value 227 the propagation SUITS. Of 
the signal. The LOPs generated from a system such as Del 
Norte Trisponder (used for this thesis) are concentric 
m jJ from the shore stations. The Del Norte system 
operates at a frequency of 9 GHz with a pulse repetition 
interval (PRI) which can be selected from 304 to 998 
microseconds. Thirty-two consecutive valid time-difference 
measurements are processed to provide each range output 
[Ingham, 1984, pp. 37-38]. 

The geometry of a range-range position fix using two 
shore stations is graphically illustrated in Figure 2.1. 
Basic trigonometric relationships are employed to determine 


the position of the survey vessel at point P using the 


T7 





Figure 2.1 Range-Range Positioning Method. 


measured ranges R, and Ry from shore control stations 1 апа. 
2, respectively. 

The computation of the coordinates of the survey 
vessel, Xp and Yp, at point P is as follows [Davis et al., 


1981, p. 876]: 


рак meses ye (2.1) 


where X], Y] are the coordinates of shore station 1l and X», 
Y, are the coordinates of shore station 2. Since all цн 
Sides of the triangle are known, angle q can be computed 
from the law of cosines equation 
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сг реке - R5-)/(2 R, D)] (959) 


Йй 3cimuth otothewlbine from station l to station 2 can be 


determined by the inverse computation equation 


41,2 - arc tan КЕТЕР 5-17)! (2.3) 


The azimuth from shore station 1 Ко the vessel P is 


mrp = 012 * a Ca 


Lf Xj and Y, are the known coordinates of the left-hand 
station when looking from the vessel toward the base line, 
then the coordinates of the vessel at P can be computed by 


the equations 


Xp = R1 sin ( 01,Р ) + X1 (2.5) 


Yp = Ry Cos ( 01,Р ) + 24 | (226) 


апа the angle of intersection of the two range LOPS can be 


computed by the equation 


ЕЕЕ > arc cos| (D> - Rugs Ro-) / (2 Ry Ro)] (2.7) 
Determining the position of a vessel uSing multiple 
ranges from three, or more, Shore stations will involve 
additional modeling since a redundancy of LOPs are 


generated. Chapter IV will address this problem in detail. 
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2. Azimuth-Azimuth Positioning Method 
An azimuth-azimuth positioning method is a direct 
extension of the standard land surveying practice of azimuth 


intersection. Ingham [1984, p. 16] states 


This [method] is usually not to be recommended, since 
complex arrangements are required to ensure that both 
angles are measured at the instant of the fix marked on 


the echo-sounder record, and to communicate the fix data 
to the boat. 





Figure 2.2 Azimuth-Azimuth Positioning Method. 
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Although true in some instances, good radio communications 
will alleviate the problem. Another option is to use a 
flagman aboard the survey vessel to coordinate the 
activities. 

The БЕ оос positioning method uses two 
theodolites, set up over known control stations, which are 
Sighted upon the survey vessel as it moves along the 
sounding line. When done correctly, the method produces a 
very accurate position at the intersection of the visual 
LOPs. The azimuth to the vessel is obtained by initializing 
on a known control station and then measuring the clockwise 
angle to the survey vessel (Figure 2.2). If theodolites 
occupy the known geodetic Shore stations 1 and 2, and 
station 0 is the station which is initialed upon, then the 
coordinates of point P, Xp, and Yp can be determined by the 


point-slope equations [Davis et нар 2711. 


Үр = Үү + (Xp - кер cot ( 0] P ) (222) 
and 
Yp = Үл + (Xp > Хо) cot ( 02 P ) (2 1395 


Subtraction of Equation 2.8 from Equation 2.9 yields 


Хр - [(Y4 - Yo) - X cot( aj p ) — 
+ X5 cot( a5 p )] / 
[cot( a5 p ) - cot( a4 p )] 
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Yp = ә ae (Xp 2 Хо) cot ( C OE ) (2,148 


Figure 2.2 shows 


Q1.p ^ 01,0 "eM (2.120 
and 
02 р ^ 02,9 * #2 (2 BB 


The use of MLOP obtained from multiple azimuths to 
determine the position of the vessel is discussed in detail 


by both Kaplan [1980] and огіуа ОИ 


a 


3. Hybrid Positioning Methods 

A hybrid positioning method (i.e., range-azimuth 
method) is a frequently used positioning method on 
near-shòre large-scale Surveys. The method provides an 
excellent geometric determination of the vessel's position 
since the two LOPs always intersect at 90°. One LOP will be 
generated from a measured electronic range to the vessel. 
The other LOP is the observed azimuth to the survey vessel. 
A concentric configuration of the ranging and angle 
measuring devices is commonly used on large-scale surveys. 
À concentric configuration is not always possible due to 
terrain and the location of the control stations. In that 


case, one of the measurements must be made from an eccentric 
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Еисссег 2.3 Range-Azimuth Positioning Method. 


Eon ТЕ (Figure 2.3). In practice station l and station 
ЈЕ are within about 2 m of each other and are occupied by a 
theodolite and a remote ranging unit, respectively. 

Knowing the coordinates of stations 1l and lE, the 
coordinates of point P may be computed by the equations 


[Wallace, 1971, pp. 46-48] 


Xp = Хүр + Ry sin( ПЈЕР ) (2.14) 


Үр = ТЕ = Ry cos ( О1Е,Р ) (424:5554 
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мһеге атр р is derived from the equation 
3 


sin( aig ,p ` @1,P ) = [1 / Ryd Uk SD) (2:15 
cos ( СР ) - ( 201 536415 Jesum 01,Р Ji 


which reduces to 


Q1E,p = @],p * arc біп (((Х1-Хір) cos( а] Р ) (2.112 


- (ҮЛҮЛ sin( С1,Р ИН R4) 


B. ERROR THEORY AND ERROR PROPAGATION 
The present NOS Hydrographic Manual [Umbach, 1976, p. 
1-4] states А ы error on any fix shall seldom 
exceed 1.5 mm at the scale of the survey. Normally 1.0 mm 
is allocated for navigational system error and 0.5 mm for 
plotting error. Thus, for a 1:5,000-scale survey the total 
allowable navigational error is 5 m. The rated accuracy of 
most positioning systems used approaches this allowable 
error. Errors must be identified and dealt with 
accordingly. 
1. Types of Errors 
Types of errors as applied to hydrographic surveying 
are described [Davis et al., 1981; Greenwalt and Schultz, 
1962; Heinzen, 1977; Kaplan, 1980; Mikhail, 1976; Perugini, 
1984 21200 
a. Blunders 
Blunders, also known as gross errors or 


mistakes, are mainly due to equipment malfunctioning or 
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observer carelessness. Blunders usually are of large 
magnitude in comparison to other errors. Unlike land 
Surveying, hydrographic surveying is a dynamic situation and 
the existence of blunders in field data is commonplace. 
When using a theodolite for position control of the vessel, 
the observer rarely has time to check the pointing to the 
moving target (survey vessel). Also the acquisition of 
multiple readings and repeating measurements independently 
to check for consistency, which are standard procedures in 
land surveys, are not feasible. On any particular fix, an 
observer has just one chance of an accurate sight to the 
vessel. If any doubt exists as to the quality of the 
Sighting, the совла Should be rejected. 

Blunders can occur when uSing electronic 
positioning systems due to malfunctioning of the equipment. 
A problem inherent with short-range systems, such as Del 
Norte Trisponder, is multipath effects (reflected wave from 
the ocean surface as well as the direct transmission) 
(Munson, 1977, p. 4]. Fading, dropout of the ranging 
Signal, and reflected signals are commonly observed during 
field work due to the omnidirectional antenna aboard the 
Survey vessel. 

Detection of blunders is one of the daily 
functions of a hydrographer. Automated data acquisition 


systems have the capability of "flagging" blunders on-line 
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through software engineering. During nonautomated 
surveying, a plot of the position fixes is maintained to 
check the accuracy of each fix as data are acquired. Most 
blunders occur during the nonautomated phase of data 
acquisition due to the increased amount of hand logging. 

Blunders can be minimized by using automatic 
electronic systems for data acquisition, such as the NOS 
HYDROPLOT Data Acquisition System (DAS) or the Racal Decca 
AUTOCARTA II system. Also, all blunders must be detected 
and resolved prior to the final processing of hydrographic 
data. Normally, repeated checking of the data will minimize 
the amount of blunders. 

b, суз ега Errors 

Systematic errors follow a definite pattern and 
are generally constant in magnitude and sign throughout a 
series of observations. The system causing the pattern may 
be dependent on the instrument, or atmospheric effects 
[Mikhail, 1976, p. 67]. In hydrographic surveying, 
systematic errors must be determined and data corrected 
prior to the final processing adjustment of all the survey 
records. 

Determination of systematic errors in 
hydrography is accomplished through the calibration of the 
surveying instruments. Calibration is the process of 


comparing the observed instrument value with a known 
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Standard. The difference between the observed and known 
values is the total systematic error present in the system. 
Once the systematic errors are determined, a "corrector" is 
applied to the observed values which yields the "true" value 
of the observations. A corrector is equal in magnitude but 
opposite in sign to the systematic error. 

À systematic error existed within the Del Norte 
Trisponder microwave ranging equipment used as the 
positioning system for data acquired for this thesis. The 
standard field procedure for this equipment was used to 
compare ranges obtained from the Trisponders to a measured 
geodetic base line. Differences were recorded and applied 
during the post processing phase of the survey as correctors 
to the ranges. Daily calibration checks were made while the 
vessel was underway to ensure that no-major changes had 
developed, since noticable calibration drifts can occur with 
time and variation in range [Munson, 1977, p. 4]. Base-line 
calibrations of short-range systems are normally performed 
before, during, and after the completion of data acquisition 
on a hydrographic survey. 

Numerous systematic errors arise when using 
theodolites for range-azimuth and azimuth-azimuth 
positioning on large-scale surveys. Proper adjustment and 
alignment of the instrument will eliminate a number of 


errors. Mikhail [1976, p. 71] lists some of the common 
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sources of systematic errors associated with the use of 
theodolites. Initial pointing and eccentricity are two 
additional sources of errors introduced in hydrographic 
surveys. 

An initial pointing must be made to a known 
geodetic station since the observer is actually measuring a 
direction to the survey vessel. Normal practice in the 
field is to adjust the plate setting to a number other than 
000°00'00", usually about 000°00'10". This value must be 
recorded in the record volume of the survey and applied as a 
corrector during processing. Should this initial corrector 
be ignored, or applied in the wrong fashion, a systematic 
error will be present in all of the observations. 

Eccentricity frequently occurs when using a 
range-azimuth positioning method as it is very difficult to 
erect a remote ranging unit (i.e., a Trisponder) and a 
theodolite directly over a horizontal control station. 
Unless two different height tripods are erected over the 
station, or a platform built to hold both instruments, an 
eccentric horizontal station must be established to locate 
one of the survey instruments. If this offset position is 
not accounted for, then a systematic error (varying in 
magnitude with the range and direction to the vessel and the 
eccentricity of the station) will result. During the data 


acquisition phase for this thesis, the ranging unit used for 
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observing the range to the vessel was erected over an 
eccentric station at three of the four horizontal control 
stations. 
с. Random Errors 
Random errors account for the remaining 
variation in measurements after all ЖЮН have been 
discovered and removed, and the observations corrected for 
the known systematic errors. Random errors cannot be 
modeled mathematically but must be modeled under the laws of 
probability. The random error is a random variable which 
can take on a number of possible values depending on the 
probability involved. 
They result from accidental and unknown combinations of 
causes beyond the control of the observer. Random 
errors are characterized by: (1) variation in sign -- 
positive and negative errors occuring with equal 
frequency, (2) small errors occurring more frequently 
than large errors, and (3) extremely large errors с 
occuring [Greenwalt and Schultz, 1962,.p. 2]. | 
Statistical treatment of random errors relating 
to hydrographic applications has been given by Waltz [1983] 
and Perugini [1984]. 
2. Accuracy of Large-Scale Survey Methods 
Two LOPs from completely automated electronic 
Systems can rarely meet the required accuracy for a 


1:5,000-scale survey due to the standard error associated 


with each LOP [Perugini, 1984]. Common survey practice has 
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been to use hybrid methods, such as range-azimuth, to reduce 
this error by introducing an accurate direction measured 
with a theodolite which will have a 90° intersection with 
the range LOP. Also completely visual methods can be 
employed, such as azimuth-azimuth or three-point sextant 
fixes, for the position determination of the survey vessel. 
All of the methods employed contain unique standard errors 
associated with the equipment utilized. This section 
discusses the rated accuracy of the equipment used for data 
acquisition for this thesis. 
a. Errors in Ranges 
The total positional accuracy "shall seldom 

exceed 1.5 mm at the scale of the survey" [Umbach, 1976, p. 
l-4]. Munson [1977, p. 2] has equated "seldom" to be 90% of 
the time which results in a 1.645-sigma value of 4.5 m on a 
1:92,000=5сате зотуеу. __ ней сре васе == 

Positioning system accuracies are most А ЫЙ, stated 

in terms of error along a line of position. Since the 

total position error will always be at least 2 

greater than this, an acceptable positioning system must 

measure a line of position to at least a lo level of 

4.5 m ЗЕ = 3.2 m for a 1:5 ООП T 

The Del Norte Trisponder system used for this 

thesis has a manufacturer stated accuracy of +1 m. Recent 
stability testing by NOS on the Model 520 Digital Distance 
Measuring Unit (DDMU) has shown errors in the order of 3 m 


over various ranges [Whitsell and Berstis, 1983, p. 4]. The 
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study also found that the range error was a function of low 
Signal strength which could be simulated by signal 
attenuation. This is a major problem because 

There is no indication in the Model 520 DDMU display of 

data output when low Signal strength conditions are 

encountered. Reduced signal attenuation will occur 

during hydrographic survey operations due to dynamic 

movement of the master antenna on the survey launch, 

inclement weather conditions, and the presence of 

multipath propagation zones. [Whitsell and Berstis, 

1983, p. 5] 

A 3-m standard error will be used for all range 
LOPs acquired during the field work for this thesis with a 
Model 540 DDMU. This unit is similar to the Model 520 DDMU , 
the only difference being that.the 540 has a four-range 
digital display versus a two-range digital display on the 
520. A fixed value over all ranges was chosen rather than 
one which varies with distance, as is normally defined for 
precise Electronic Distance Measuring (EDM) equipment in 
land surveying, to simulate common field practice. AS a 
result, the error of each LOP at a short range (1,000 m) 
will have a higher bias than a LOP at a longer range (10,000 
m) due to the fixed error. 
b. Errors in Azimuths 
One of the objectives of the work done by Waltz 

[1983] was to determine the pointing error of the Wild T-2 


theodolite when used for range-azimuth hydrographic control. 


His rigorous treatment of the subject found that an 
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estimated linear error for the pointing of a T-2 was 1.3 п. 
Actually the pointing error is a function of the angular 
resolution of the instrument plus the ability of the 
observer to track the sounding vessel. Both are a function 
of the distance to the vessel. For azimuth-azimuth 
positioning methods, the errors in positioning actually 
depend upon the distance, the angular resolution, and the 
angle of intersection of the LOPs [Heinzen, 1977, p. 55]. 

An important conclusion made by Waltz (|1983, р. 
81] was 

There exists, for any angle measured with a T-2, a time 
lag of about one second between angle observations and 
any measurement made aboard the vessel, including both 
automatic and manually recorded depth and range data. 
There is then an associated position error for these 
measurements, the magnitude of which depends upon vessel 
speed, which was about two meters for the four knot 
Speed used in this experiment. 

With this fact in mind, a standard error 
associated with azimuth LOPs was computed for the data 
acquired for this thesis. Hand plots of the position fixes 
show that the speed of the survey vessel was approximately 5 
knots (2.6 m/s). Thus, the l-second time lag would amount 
to as much as an additional 2.6 m of error in position due 
to the vessel velocity relative to the shore station. 
Consider an example which shows positional error as a 
function of the range and the angular resolution of the 
instrument (Figure 2.4). Using the standard conversion 


formula 


32 





Figure 2.4 Conversion of Angular Measure to Arc Length. 


р = Sp | R pos 


where b is the arc swept by the vessel (2.6 m), R is the 
range to the vessel in meters, q is the angular error in 
seconds, and p Ea conversion factor (206265"/rad). If R 
Ээ then a S (206265"”)(2.6 m. / 1,000 m) = 536" = 
ИСЕТ СЕ, 16 К іс 2.000 ш, сђеп а_= 4.5': if R is 
3,000 m, then q = 3.0'. 

A total standard error for a LOP measured by an 
azimuth for this thesis was assumed to be 4m. This value 
takes into account a 1.3-m pointing error plus the 2.6-m 
error due to the l-second time lag at a vessel speed of 


approximately 5 knots. This was done to duplicate standard 
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processing techniques for azimuth data. By assuming a fixed 
error over all ranges a reverse situation holds for LOPs 
from azimuth measurements when compared to the similar 
situation for range measurements. At longer distances the 
measured azimuth will be more accurate than at shorter 
distances because the aspect of the vessel does not change 
as rapidly. By using a fixed error in position, the 
azimuths obtained at close .range will have an underestimated 
error whereas those observed at longer ranges will have an 
error which is overestimated. 
3. Error 2004-0881 

Since random errors will exist in the field 
measurements, then the errors in the computed positions of 
the survey vessel can be obtained by the law of propagation 
of errors. Error propagation is the evaluation of errors in 
the computed positions as functions of the errors in the 
field шоогоо Ав 21 previously, random errors 
follow statistical laws and can be dealt with through 
probability distributions of these errors. 


Let the function Y be defined by 


Y = £(X1,X5,X4, . ° чай (2 Жар 


np 


in which Rp. DOMENUS Ke uncorrelated 
measurements. Thus, the law of propagation of errors 


[Mikhail and Gracie, 1981, pp. 148-178] can be shown as 
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oy^ - (OY /| 0X )* ox,* « ( ӘҮ / 2%; )4 су, (2.20) 


2 
p ( Qu / Qu йн ox ^ 


Consider the problem of defining the position of a 
survey vessel at sea. The goal of the positioning methods 
previously discussed is to accurately determine the computed 
coordinates (Xp. Yp) of the survey vessel. Each of the two 
coordinates are independently a function of the measurements 
used to compute the coordinates. For example, when using a 
range-range positioning system the coordinates are a 
function of the measured ranges (Ry and Rj) and the known 


station coordinates on shore such that, 


Xp = Ғ1(Е1.,К2,Х1.Х2) ; (2.21) 


where КШ ҮҮ and X5, Y are the known coordinates of the 
shore stations l and 2, respectively. The law of 


propagation of errors takes the following form 
gue Q ORw/OR. )? са.“ +  ( OXp/üR; )? св, (0,22) 
* ( 0Xp/0Xi )* oy,^ * ( OXp/OX) )* оху 


су“ = ( OYp/GRi 0^ og, ^ *  ( OYp/OR, )* св, (2.24) 


+ ( дҮр/дҮ1 )2 бү,“ AN QYp/ OY» y су,“ 
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The computed position (Xp, Yp) from a range-range 
positioning system was given by Equations 2.1 through 2./7. 
Partial differentiation of the terms in Equations 2 29 420 
2.24 yields Equations 2.25 through 2.28 where 


then 
ОХр / 081 - біп( 0] P 2. [cos a1 P )] (2.2) 
( 1 - u2 ) ~ (ке - ~ a 
дХр / OR 7 [cos( a4 p )] ( 1 - u? ) V? (2.26) 
(Rə / D) 
lp / ORj * cos( aj p ) * [sin( aj p )] (2.20 
(1-4«2)7М2 ((в12» в,4- 22) / (2 кү) 
дур /| OR 7 [sin( a4 p )] ( 1 - v? ) !/? (2.28) 
баса 2 | 


Since the X and Y coordinates of the shore stations are 


constant then 
дХр/дху = дХр/дх; = дур/дуј = дур/ду; = 0 (2.998 


Previously, it was shown that the standard error of 
the Del Norte Trisponder system, Ор, has been determined to 


be 3 m, regardless of range, so 


ORY = бұ, = бр = 3 m (2.30) 
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Then, the standard error of each of the computed coordinates 


will be 


бу = ср / (ОЖр / OR,)* (дїр / 285)2 (2.31) 


and 


Buc ср 200055) 25172 > (гүр / око)? (2.32) 


For range-azimuth positioning systems, the 
coordinates of the vessel are a function of the measured 
range and azimuth on each fix plus the X, ¥ coordinates of 


the range and azimuth station 


and 
Yp = 55 (К1.01.Р,Ү1Е»Ү1) (254) 


so that the equation of the propagation of errors takes the 


form 


Sx^ * ( OXp/OR, )^ og,^ * ( OXp/oui,p 0^ og, ,? * (2.35) 


^ 2 2 2 
( OXp/oXig )° оу,„© * ( дХр/ду, )° оу,“ 
oy? S 0Yp/ 0R1 )* OR, + ( дҮр/дал р )2 с а (2 6) 


2 


^ 2 2 2 
( дҮр/дҮүр ) бүүр ин OYp/ Oy. ) Ou 
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The position coordinates (Xp Yp) were given in Equations 
2.14 through 2.17. The value ОЕ отр p (Equation m 


repeated here (Equation 2.37) for continuity. 


ü]g,p 7 arc sin ([(Х]-ХјЕ) cos ( а] р ) (2.37) 


- (Үү-Үүр) sin( а] р мг @Т,Р 


Partial differentiation of the terms in Equations 2.35 апа 


2.36 where 


u = (1/R,) [(Ху - Хур) соз(а) р) - (Ү) - Те) sina, p)] 


yields 
ОХр / 081 = 510(01Е р) - u | со8(атр ур)! (2.39 


(1 - и2 )-1/2 


дХр / дај „ = (сов(ајр р)Ј (Ер - ( 1 - u? )y-1⁄2 — (2.39) 


(Xy - Xig) [sin(aj,p)] * (Y, - Ті) соз(ај p)} 
дҮр / ОК1 - coS(aig p) оп [sin(ayp p)] (2.40) 


о а. 


Yp / дау ь = Їзїп(аукЕ,р)] {[( 1-ч°^ )7+/5 - ву] (2.41) 


OG š Үүр) |соз(01 р)) T (Хү я Хүр) 51п(01 р)! 
Since the X and Y coordinates of the range and azimuth shore 


Stations are constants then 


QXp/OX1r = 2Хр/0Хі = дҮр/дҮүр - дҮр/0Үү - 0 (2.42) 
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Substitution into the propagation of error equations 


(Equations 2.35 and 2.36) results in 


суг = ( e ORI уе сву + ( дХр/дај р ЈЕ O41. P^ (T) 


сүг = ( 0Yp/OR, 77 св + ( дур/дај,р 0 бил РО (2:44) 


and the standard error Ox and Oy of the X and Y coordinates 
is obtained from Equations 2.43 and 2.44. 

The standard error of position, Ор» was used as a 
reference to classify the accuracy of the data acquired for 
this thesis. Rather than classify the accuracy using the 
‘standard error ellipse or confidence circles, the бр value 
was chosen due to its ease of computation. The value of the 
Шанаға position error [Saxena, 1972, p. 15] is defined as 


the square root of the sum of the squares of the propagation 


variances of the system. Mathematically 


а 2 2 
Op = oy Ov (2.45) 
where ane and су“ are computed as shown previously for 
range-range LOPS and range-azimuth LOPs. The standard error 


ОИ О Ен от са in em, for azimuth LOPs (Mueller, 1979, p. 


621] 15 


og = ( 0.485 / sin B ) o, Бе + Do" (2.46) 
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where B is the angle of intersection of tne two БОР ТН 
the distance from point 1 to point P in km; Dy is the 
distance from point 2 to point P in km; and Bi is the 
angular error in seconds of arc. 

For the overdetermined case of MLOP, the g, value 


P 
(as derived by Equations 4.28, 4.29, and 4.30) will be 


Ото Охх (2.47) 
Cae Or Qyy (2.48) 
and 


Ор ^ J Ox2 * oy2 7 og J Qxx * Qyy (2252 


where QE is the standard error of an observation of unit 
weight, Qxx and Qyy are the diagonal elements of the 
variance-covariance matrix. 

Thus a tag can be placed on every position fix 
acquired саа whatever system, or combination of systems, 
used to determine the coordinates of the survey vessel. 
Using this standard position accuracy, a comparison between 
fixes determined by different systems and methods can be 
compared so that the most precise combination of ranges, 


azimuths, or hybrid methods can be analyzed. 
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111. DATA ACQUISITION PROCEDURES AND SYSTEM PERFORMANCE 


The field work required for this thesis involved 
conducting a hydrographic survey while using MLOP for 
sounding line position fixing. Two types of positioning 
data, ranges and azimuths, were acquired simultaneously from 
four shore control sites. The data were then analyzed using 
various combinations of ranges and azimuths to determine a 


more accurate positioning system. 


A. DATA ACQUISITION PROCEDURES 

Е: диско in the NPS Hydrographic Sciences curriculum 
conducted a basic hydrographic survey of southern Monterey. . 
Bay as part of a field experience requirement. The survey 
was accomplished using procedures similar to those used by 
NOS for nearshore surveying. A chartered 36-foot, 
twin-engine Uniflite е. was the platform used for the 
Survey. During the week of 26 November 1984 to 30 November: 
1984, the boat was positioned using four ranges and four 
azimuths from known third-order horizontal geodetic control 
stations on shore for each position fix of the vessel. The 
ranges to the vessel were determined from four Del Norte 
Irisponders and were automatically recorded by a Racal-Decca 
AUTOCARTA II data acquisition system. The four azimuths 


were observed with Wild T-2 theodolites. Data acquisition 
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equipment was supplied by Racal-Decca, NOAA, and NPS 


(Table I). 


TABLE I 


EGULEMENT LIST 


Aboard M/V Silver Prince 


Model 540 DDMU S/N 211 (Racal) 
Master Trisponder S/N 3014 (NOAA) 





Onshore Control Stations 
STATION REMOTE TRISPONDER T-2 (S/N) 


SQUARE 1984 Code 72 14405 (NPS) 
S/N 2819 (Racal) 


CONK 1984  . Code 74 | 51642  (NPS) 
S/N 2822 (Racal) | 


USE МОМ 1978 Соде 76 30504 (NPS) 
S/N 3004 (NOAA) 


GEOCEIVER 1982 Code 78 14482 (NPS) 
S/N 2986 (NOAA) 


Geodetic control for the survey consisted of third-order 
monumented stations and established eccentrics (Table II). 


The eccentric stations were within 2 m of the main station 
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and marked by masonry nails driven into the concrete. The 
geodetic position of each eccentric station was computed on 
a HP 9815 calculator using National Geodetic Survey (NGS ) 
geodetic programs. Initial pointings for the T-2 
theodolites were to MUSSEL from stations SQUARE, CONK, and 
USE MON; and to MONTEREY RADIO STATION KMBY MAST from 


station GEOCEIVER. 


TABLE II 


GEODETIC CONTROL POSITIONS 


Station X (m) (m) Latitude Longitude 


SQUARE 7974.86 .43 3637 "07.175" 121951 0072 4 
SQUARE ecc 7976.10 ,89 36%37"07.140" 12195177007 Ж 


CONK 6978.19 .77 36936 32.130" 1219578 ИВ 
СОМК есс 6976.70 .53 36936"32.100" 121%51 40,4 БИШ 


05Е МОМ 5598.98 .76 36736'04.685" 121"52 #35 ШИН 


GEOCEIVER 4371.50 .28 36936"32.512" 121953 254 ШЕ 
СЕОСҮК есс 4372.99 .97 36%36"32.5107%7121953 ЖА 


MUSSEL 52208 ,23 36937118.151" 121954711 ШИШ 
KMBY MAST 3641.23 .23 36°36'56.789" 121?53 549 ED 





The data acquisition procedures were similar to normal 


survey operations for a large-scale hydrographic survey 
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Memes visual control with the addition that on each fix four 
ranges were recorded automatically by the AUTOCARTA II 
system. The four theodolite observers trained their 
instruments on the vessel as it moved along the track line. 
Position fix marks were communicated over voice radio and 
the T-2 directions to the vessel were manually recorded by 
the observers. 

To eliminate systematic errors associated with the Del 
Norte system, the Trisponder units were calibrated over a 
measured geodetic base line 2565.897 m long (CONK 1984 to 
REY 1984) on 24 and 25 November 1984. A base-line 
calibration prior to survey "— is a standard field 
procedure for this type of equipment. Daily dynamic system 
checks were made EH S NOS procedures [Holder, 1983]. 
The она "drifts" in the Trisponder units were recorded 
from an ending base-line calibration on 7 December 1984 and 
were applied as range correctors during the data processing 
phase of the experiment. 

Weather conditions during the 2 days in which horizontal 
directions to the vessel were recorded was fair. The 
theodolite observers had difficulty in maintaining an 
adequate initial pointing due to winds gusting to 30 knots 
which caused vibration of the tripods. Initial checks were 
made at the end of each track line and recorded. Observed 


changes in the initial directions were recorded and applied 
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during the data processing phase to eliminate any systematic 
errors in the angles. On the first 2 days of field 
operations, a total of 1/6 positions, each with four ranges 
and four azimuths, were recorded. An additional 175 
positions with four ranges were acquired using the AUTOCARTA 


II system during the last 3 days of the Ese dc USE 


B. DATA ACQUISITION SYSTEM 

The AUTOCARTA II (Figure 3.2) is a portable automated 
hydrographic data acquisition system. Desired vessel tracks 
can be preprogrammed into the computer to alleviate the 
necessity of steering the vessel along range arcs for track 
control. A real-time graphic copy of the vessel's actual 
course was plotted from Autocarta II output on a Houston 
Instruments Model DP-3 Plotter. A Model 540 DDMU recorded 
the four ranges at ТЕ intervals along the trad line. 
Input and output: comnunication to the AUTOCARTA II were via 


a modified TI-743KSR keyboard/printer. 


C. PERFORMANCE OF DATA ACQUISITION SYSTEM 

The operational аже experienced during the data 
acquisition phase of the field work could be described as 
"standard" for a hydrographic survey of this type. To 
ensure adequate backup was available for Trisponder units, a 
total of seven remote and two master units were on hand. 


Two of the remote units did not responded when interrogated 
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Data Acquisition System. 


Figure 3.2 


during the base-line calibration of the system and had to be 
removed for the entire study. New deep-cycle marine 
batteries were used as power to the 24-volt Trisponder 
system. These batteries were charged each night to ensure 
sufficient power during the survey. | 

A cable linking the AUTOCARTA II to the keyboard/printer 
was not shipped with the equipment. Therefore, the data 
were acquired in the normal non-automated fashion by hand 
recording the range rates as they were displayed on the DDMU 
on the first day of the survey. The needed cable was 
available by the next day and the complete automated system 
was operational for the rest of the survey. 

The AUTOCARTA II was set up to interrogate the remote 
units on shore from left to right when facing shoreward from 
the vessel which is standard NOS convention. On the first 
day of data acquisition, day 332, very erratic range rates 
were being observed on the DDMU from station USE MON 1978 
(remote code #76). The update rate was set at l-second 
intervals and over half of the time either a totally 
ambiguous rate or no rate at all was recorded on the DDMU. 
In normal field operations, a remedy to the situation is to 
raise or lower the remote unit to eliminate the possibility 
of "null zones" and "skip zones." The remote unit on 
Station USE MON 1978 was raised and lowered with no change 


in the performance of the Trisponder observed. 
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IV. MATHEMATICAL MODEL 


A. LEAST SQUARES MODEL 

To determine the "best" estimate of a position, 
redundant observations are necessary so that the method of 
least squares can be used. The least squares adjustment by 
observation equations (variation of coordinates) was 
selected for this thesis (rather than least squares 
adjustment by condition кесі тий for reasons of simplicity 
ШІ сіатісу [Bomford, 1980, pp. 127-128]. 

A basic discussion of the method of least squares by 
observation equations is given below [Bomford, 1980; Mikhail 
and Gracie, 1981; Cross, 1981; Heinzen, 1977; Бог ас, 1300: 
Saxena, 1972; Silva, 1982}. 

l. Unweighted Equations 
The equation derivations will be initially 
accomplished for measurements of equal accuracy. 
Ја represent independent аена quantities; vj; 
ТВ наас of the observed quantities; and x, y, and z be 
the unknown coordinates of a point in a rectangular 
coordinate system. Then for each observation, i, an 


observation equation is formed such that 


Ша + Vi = Ши) (4.1) 
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where f; is a linear or non-linear function. When i 15 а 
non-linear function, a Taylor series is formed about the 


and z. such that the second 


approximate values of xg, yg O 
3 


and higher order terms are small enough to be neglected. 


The observation equations then become 


І. Ру; SUES CRS Т ат” 222. 
Og 
С е шы ы (47752) 
+ a,dx + b;dy + cidz + 
EE ДА = I Е ИТ 97А which is the difference between 


the observed quantities and the computed approximate values 


of the unknowns, then 


Vi = a; dx + b; dy + c; dz Pire foin 24 (4.4) 
where 

х = хо + Ях (42458) 

y NE SONGS; (4.6) 

22-02 clo (4.7) 
and 

a, = oE sp s ш> (&.8) 
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so that in matrix form the observation equation 4.4 becomes 


ЕЕ А)- L (4.9) 


which ultimately yields the normal equations 


NX - U = 0 | (2210) 
where 

N = АТА | (4.11) 
and 

U = ATL | (4.12) 


Then substituting into Equation 4.10 and rearranging terms 


gives 
x = N`lU = (ATA) 1 ATL (4.13) 


which is valid for pices нене of equal weight. 

In actual surveying, computed observations differ 
from actual observations by a value known as a residual, v. 
The residual is the difference between an observed value, 


аза ВЕ estimate ofthe true value, x, so that 
5 


X = x; + У. (4.14) 


The basic condition of the least squares method is 
that the sum of the squares of the residuals equal a 


minimum. Therefore, 


SUL 


Y 2208 = minimum (4.15) 


where 
> 227 - Б + E и ас 
for values of equal weight. 
2. Weighted Equations 
When observations are made utilizing various methods 
and systems, weights must be applied to account for the 
lower standard error, o, and lower variance, | of the 
more precise observations. The relationship between weights 


and с? is 


w = k / o ` (4.10 


where k is a constant of proportionality having the same 
value for all го in one system when the different 
standard errors of the observations are uncorrelated and of 
the same precision, as is the case for range-range, 


azimuth-azimuth, and range-azimuth positioning systems. 


Thus 
Then 
21 2 4.18 
W; T (ёс (4.18) 


D 


meh іп matrix form is the diagonal matrix 
1/62 
1/g5* 


М = | (4.19) 
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For weighted observations, Equation 4.15 becomes 


£ 202 = minimum (4.20) 
where 
> 2707) = w1v1* + ud гэ: Wave 


ӘП in matrixəform 
viwy = minimum (4.21) 


which after substitution and differentiation with respect to 


X becomes 
x = (АТНА) 1 АТ, (4.22) 


This system of equations is then solved for the 
matrix X elements which are then added to a provisional set 
of coordinates at a point, yielding a new pair of 
coordinates. This process is iterated until the values of 


the X elements become small enough to meet a Specific 


2 


tolerance. At this point the method is said to have 
converged. What remains is the "best" estimate of x and v, 


A A 
designated as x and у. 


B. PRECISION OF ADJUSTED VALUES 
The first right-hand term of Equation 4.22 is the 


variance-covariance matrix, 
(АТА) У = 1 = 0 (4.23) 
and it can be shown that, 
411 412 
ЗӘЛГ 222 


Q = (ATWA) E . QUE 


dnn 


The variance-covariance matrix of the least squares estimate 


О (3), 15 


с($) = 02 (АТНА) = ои (a 
4 = 
Ox Oxy 
2 
бух бу 
g(x) = 
2 
ба 
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Here тоо lS known as the variance of ап observation of unit 


weight represented as 
б^ - viwy / (n 2 m) 210) 


The standard error of an observation of unit weight, gg 15 
5 


therefore 


со = ./ Viw / (n - m) (4.27) 


where n is the number:/of observations, and m is the number 


of unknowns. The term (n - m) is known as the degree of 


freedom. 


Consequently, the standard error of the adjusted 


T 


coordinates, Gs and g are given by; 


e 
ox = oo / qx - (4.28) 


Ov = Об ~ Чуу (0229) 


апа the standard position error is 


Ор E. сұ” " ore к с A О + Оуу СОЗО) 
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Finally, the semi-major (oz), semi minor (op): qmi 


orientation (0) axes of the error ellipse are given as 


Sa = Oo V Чпах (4. 


Ob = Со y Qmin | (4. 


_ where 
Qnax,min = Emm (Qus i 0уу) (4. 
+ 4 2 2 
£ (1/2) J(Qxx-Qyy)^ + 4 Оху 
and 
0 = (1/2) arc tan [ 2 Qy / (Qu, - Оуу) ] (4. 
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V. DATA ANALYSIS AND RESULTS 


A. DATA PROCESSING PROCEDURES 

Data acquired in the field were manually logged into 
data files on the IBM 370/3033AP mainframe computer at NPS 
via the VM terminals. These data were processed using 
computer programs (Appendix A) written to be compatible with 
the VS FORTRAN compiler and used many of the existing 
library subroutines. The main program which computes the 
least Squares position from three- or four-range LOPs, 
SLEAST (Appendix B), was composed by Lt. Nick Perugini, 
NOAA, and the author. The article by Cross [1981] covers 
the subject in detail. The least Squares programs for the 
azimuth and hybrid MLOP methods are FORTRAN programs LSQAZ4, 
R3*AZ, R2*AZ2, and RR+AZ (кррецат # A). These programs are 
modifications of programs SILVAl and SILVA3, written by 
БЧБ S [1982]. Computation of coordinates and accuracies 
from two LOPs were made by modifying the program UCOMPS 
(Appendix A) into RR2LOP, AZ2LOP, RAZ2LOP. 

Processing procedures commenced by correcting the 
observed ranges and azimuths for previously determined 
systematic errors prior to logging the data into manageable 
files. The least squares computer programs output the 


position of the vessel (in X and Y coordinates based on a 
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Modified Transverse Mercator (MTM) projectionl); the final 
residuals; and the computed accuracy, Op» of each position 
fix. The output from the processing programs for two LOPs 
were the position coordinates of the vessel, the angle of 
intersection of the LOPS, В, апа the запао вани и 
Positron, Ор» 45 determined in Chapter II for each position 
method. 

The geodetic positions of the shore stations used for 


hydrographic control were converted to MTM coordinates via 


the utility program UCOMPS (Appendix A). 


B. DATA ANALYSIS 

The main objective of the thesis was to investigate 
whether positions computed from MLOP were more accurate than 
positions determined by only two LOPs. Models and 
algorithms were developed so field data could be 
appropriately analyzed. For each positional fix, eight LOPs 
(four ranges and four azimuths) were measured 
Simultaneously. The eight LOPS were processed separately 


into ten different sets of positioning methods. Depending 


lrhe МТМ projection is a Universal Transverse Mercator 
projection except the Central Meridian is near the Survey 
area instead of being at a preselected meridian [Wallace, 
1971, рр. 3-6|. А Сепега1 Мегійзап оҰҒ longitude 121°53 ии 
was chosen for this thesis. The X coordinate assigned to 
the Central Meridian, the FEST, was chosen as 5,000 m. The 
Controlling Latitude, CLAT, or the distance in meters from 


the equator to some reference latitude, was computed to be 
4,050,000 m. 
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on the method, some of the ten sets were separated into 
various subsets by using different combinations of 
observations from the four shore stations. 

Let the four shore stations from which the LOPs 
Originated be designated A, B, C, and D. If four-range (or 
four-azimuth) LOPs are to be used for positioning, there is 
only one unique combination of four-range (or four-azimuth) 
LOPs which can solve the problem. Using only three LOPs, 
there exists four possible combinations of range (or 
azimuth) LOPs--ABC, ABD, ACD, and BCD. | 

The problem becomes more confusing with hybrid 
combinations of range LOPs and azimuth LOPs. A practical 
limitation was imposed by deciding that each azimuth LOP 
must originate from the same station (or eccentric station) 
as the range LOP. Thus, if three LOPS were to be determined 
from two ranges and one azimuth, the azimuth observation 
ЕЕЕ Бе made at one of the stations (or its eccentric 
Station) where the ranging unit is located. For example, if 
B issonders are set on stations A and B (or the ОООО 
stations АЕ апа ВЕ). The azimuth LOP must originate from 
either station À or station B. There are a total of 12 LOP 
combinations which must be considered when uSing three 
ranges and one azimuth, or two ranges and one azimuth 


(Table III). 
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TABLE III 


LOP COMBINATIONS FOR TWO HYBRID METHODS 


3 Ranges & 1 Azimuth 2 Ranges & 1 Azimuth 





3 Ranges Azimuth 2 Ranges Azimuth 


g OQ = g O > Ü Z > OQ Z > 
g Q g = Q = Ü > G) > Z > 





All possible combinations of the 10 data sets were 


processed resulting in a total of 56 positions being 
computed for most of the fixes (Table IV). 

Gaps exist between consecutive fix numbers in the output 
files (Appendices B, C, and F) due to spurious ranges or 
unresolved blunders which were rejected. Data were rejected 


if the standard error of position exceeded 3с for that LOP. 
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TABLE ТУ 


FILE COMBINATIONS 


Method Possible Combinations 





The following assumptions were made concerning the 
overall accuracy determination of each position: 
(a) All LOPs were independant and uncorrelated. 
(b) All data were free of blunders and systematic errors. 


(c) The standard error of range LOPs was 3 m at all 
ranges. 


(d) The standard error of azimuth LOPs was 4 п. 
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(e) The measured ranges were referred to the same point 
of time during interrogation by the DDMU time-deskew 
function [Anonymous 0 8. 


(Е) Measurements by the azimuth observers were all of 
equal weight. 


1. Range Measurement Methods 

An analysis was Ае of the acquired range data by 
the number of LOPs. The mean value of бр from each 
positioning method is illustrated in Table V. 

а. МЬОР 

FORTRAN program SLEAST (Appendix B) processed 

the three- and four-range LOPs (Appendix C) using the least 
squares method. The mean Ср value for each day of data and 


the minimum and maximum g, were computed (Table V). For 


P 
the 3 days in which range measurements were acquired, an 
average o, of 1.3 m was attained with the four LOPs and 1.5 
m with the three LOPs. The geometric accuracy of the MLOP 
fixes exceeded the standard error af the Del Norte system by 
at least 1.5 m. The slight variation in the mean бр 5 БОБ 


the three LOPs was due to the geometrical variation of the 


ranging net from each shore station. 


b. Two LOPs 
The two-range LOPs were processed by program 
RR2LOP (Appendix A) to compute the coordinates of the vessel 
plus the standard error of position, Op» for each fix. 


Analysis of the data shows a direct relationship between the 
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TABLE V 


POSITIONAL ACCURACIES FOR RANGE LOPS 





Op Op Op 

LOPs DAY MIN MAX MEAN 
(m ) (m) (m) 

& R ЭЭ? 0.3 ZG 1e. 
GR 333 0.4 дж 9 Кол 
& R 334 0.5 Du d 
в-к 332 0.1 ЗөвӨ ls 
0.04 45 1.7 

S R 333 0.04 5 15:53 
0.04 Je 29 

3 R 334 0.1 2.4 0.8 
0.4 4.5 2212 

2-8 332 4.8 7 6.0 
7.7 23567 2M. 

28 333 4.4 ga 5 
дас? 12-3 


* NOTE: For positions which could be 
than one combination of LOPs (e.g., 4 


ANGLE OF 
INTERSECT ION 


(2 LOPs) 


147? 
2-9 


28° 
Де 


computed by more 
combinations of 


3-range LOP fixes exist when 4 ranges have been 
measured), the upper line for each day shows the 
combination of LOPs giving the minimum mean gp's 


and the lower line for each day shows 
of LOPs giving the maximum mean gp's. 
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the combination 


angle of intersection of the LOPs and ср, аз would be 
expected. A wide variation in the computed Ор 5 exists due 
to the varying angle of intersections. The pre-selected 
sounding line orientation was programmed into the AUTOCARTA 
II prior to the acquisition of data. Since the purpose of 
the field work was to determine the position of the vessel 
using MLOP, the angle of intersection of each individual 
pair of LOPs was ignored. Thus at fix number 19 the 
computed Op was 237 m due to an angle of intersection of 
1.69. During normal survey operations this situation would 
not be accepted and alternate plans would be made to obtain 
coverage in that area. Herein lies a benefit of using MLOP 
Since the position of the vessel can be determined from 
selected shore stations ши yield the best geometric 
configuration. As could be expected, the required accuracy 
of at least 5 m was not met when using two range LOPs with a 
standard error of 3 m (Table V). | 

Overall the capability of acquiring MLOP 
ачк EE ME from measured ranges has proven to be a very 
viable system. Mean accuracies from 1.2 to 1.6 m were 
observed depending on the number and geometric combinations 
utilized. With these accuracies all nearshore positioning 
requirements can be met. 

Using range MLOP for hydrographic positioning 


has many advantages. The capability of running a fully 
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automated survey with systems such as the AUTOCARTA ITI not 
only increases the productivity but also reduces blunders, 
thereby increasing the quality of the survey. Once remote 
transponders are erected over the control stations the need 
for attendant operators is negated thus reducing the number 
of personnel required to acquire data. The portability of 
the system (except possibly for the batteries) is an 
additional advantage in rugged terrain. 
2. Azimuth Measurement Methods 

The use of an angle measuring instrument, such as a 
T-2 theodolite, for positional control of sounding lines is 
a laborious task. Continuous tracking of the vessel is 
difficult because the horizontal tangent screw has a finite 
range. Observers can be tracking the vessel with the 
tangent screw only to find the end of the drive mechanism 
seconds before the position fix. Thus the problem of 
tracking a vessel who's aspect is changing rapidly with 
respect to a nearby observer is magnified. Older T-2 
theodolites have an inverted image of the target which may 
confuse an observer since the vessel seems to be moving in 
the opposite direction. Also the reading of the observed 
direction is through a Secondary Scope. Therefore, an 
observer must stop tracking the vessel, align the horizontal 
plates, read the direction and record the value. The lack 


of digital telemetering capabilities requires the observed 


65 


value to be relayed to the survey vessel by voice radio. 
Obviously, there exists a potential for many sources of 
errors if position fixes are required at frequent 
(30-second) intervals. 

Using nonautomated systems (i.e., visual azimuths) 
restricts position fixing of the vessel to well spaced 
intervals (minimum of 30 seconds). Intermediate positions 
along the sounding line must be interpolated via 
dead-reckoning. In most cases, the actual track of the 
Survey vessel does not exactly correspond to the plotted 
survey line. On a 1:5,000-scale survey, this additional 
error could exceed the specified accuracy requirements. 

Although the above limitations seem to make em 
method too complicated and error prone for hydrographic 
control, it is still one of the most exact positioning 
methods available for large-scale surveys. Experienced 
observers have little dU COT in acquiring accurate data. 

a.  MLOP 

Even though a 4-m standard error was applied to 
all azimuth LOPs, the least squares position from three or 
four azimuths proved to be the most accurate method 
investigated. Examination of the tabulated бр 5 (Table 
VI) shows a mean value of 0.7 m on day 332 and 1.0 m on day 
333 when the four LOPs were used. The mean g,'s from three 


P 


LOPs ranged from 0.7 т to 1.3 m over all combinations from 
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TABLE VI 


POSITIONAL ACCURACIES FOR AZIMUTH LOPS 





Op Op Op ANGLE OF 
LOPs DAY. MIN MAX MEAN INTERSECTION 
(m) (m) (m) (2 LOPs) 
4 de 2217: 1.7 RT 
-— = 27 2.9 1.0 
D 332 0.00 2.5 0.7: * 
0.01 9.5 1.2 х 
DE ЖӘ 2055 450 0.8 : 
ONCE E 1.3 * 
2E 7 8.9 823 41° 
бозие є 216 7 19 
2 AZ тг”: — 26° 
5.7 84.2 8.6 до 


alo 


* NOTE: For positions which could be computed by more 
than one combination of LOPs (e.g., 4 combinations of 
3-azimuth LOP fixes exist when 4 azimuths have been 
observed), the upper line for each day shows the 
combination of LOPs giving the minimum теап бр 5 

and the lower line for each day shows the combination 
of LOPs giving the maximum mean gp's. 
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the 2 days of data. The accuracy capabilities of the method 
are evident. 
b. Two LOPs 

Although the least squares adjusted positions of 
the azimuth MLOP proved to be extremely accurate, the 
two-azimuth-LOPs condition suffers from the same degradation 
in accuracy as the two-range-LOPs case. The error in 
position, os in cm, was determined from Equation 2.46. 
Since the standard error of an azimuth LOP has been defined 
as 4 m throughout this thesis, ga must be evaluated. From 


the standard conversion formula (Equation 2.18) then 


с = ( b / R ) 206265 s/rad (58 


where R is the distance in m, and b is the standard error of 
4 m. Since the azimuth LOPs result in Xp and Yp, the 
coordinates of the vessel, then the distance from the 


computed coordinate to each shore control station is 


J Gq-Xp)? * (Y4-Yp)* (5 28 


d1 


|| 


апа 


42 = |/ (Xy-Xp)* + (¥9-Yp)* (528 


Rather than compute separate error values for each distance 


an average of the distance from both stations, d was used 


а” 


where 


а= ааг) о (5.4) 
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and Equation 5.1 becomes 


C oim ( 4 / а„ ) 206265 (21:22) 
which can be substituted into Equation 2.46 to compute og 
(in cm). Consequently, the computed бр (o, in cm times 

100) is a function of both the range to the vessel and the 
angle of intersection of the two LOPs. 

Analysis of all of the Op 5 attained in this 
fashion shows a mean value of about 6.3 m for all data 
acquired (Table VI). Although this is rather high for an 
error in position from two accurate azimuths, the cause сап 
be easily explained. The initial assumption that T 
azimuth LOP was In erm by 4 m at all ranges is rather 
excessive for theodolite measurements. A more realistic 
investigation of the angular error of dynamic azimuth . 
measurements should be considered in the future. Standard 
field procedure for azimuth hydrographic control is to 
ignore the stated random error EERE theodolite because the 
measured directions are only read to the nearest minute of 
arc. The only errors associated with a position determined 
by azimuths are the initial correctors and any blunders 
discovered during data acquisition and post processing. 

3. Hybrid Methods 
Numerous combinations of control configurations 
exist when ranges and azimuths are mixed together. The data 
acquired were processed by FORTRAN programs R3+AZ, R2+AZ2, 
and RR+AZ (Appendix À). 
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a. Three-Range and One-Azimuth LOPs 
Three-range and one-azimuth LOPs exhibit larger 


mean values of "5 than in the prior MLOP methods (Table 


Ор 
VII). One explanation for this difference is the geometry 
of the fix has a greater effect on the hybrid methods than 
on the other methods of position control. In тө. 
situation, one intersection is at 90° and up to four 
intersections are at arbitrary angles depending on the 
position of the vessel relative to the shore stations. If 
the angle of intersection of the range LOPs 15 too large 
(greater than 150°) or too small (less than 30°) then the 
quality of the fix may be degraded. The 11.8-m and 12.5-m 
maximum values for the Op S observed on both days (days 332 
and 333, Table VII) were found using configuration ABDD 
(three ranges from stations А, B, and D and the аз ии from 
station D). In both instances, the angle of intersection 
between the LOPs enon stations A and B was greater than 
150°. Therefore, the hybrid method is very sensitive to the 
angle of intersection of each LOP pair where a range and an 
azimuth are not taken from the same station. 
b. Two-Range and Two-Azimuth LOPs 

Analysis of the results obtained from processing 
the two-range and two-azimuth LOPs found not as much 
variability in the range of the o,,'s as in the Ор 5 of 


P 


three-range and one-azimuth LOPs. An overall mean Ор of 
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TABLE VII 


POSITIONAL ACCURACIES FOR HYBRID LOPS 





Op Op Op ANGLE OF 
LOPs DAY MIN MAX MEAN INTERSECTION 

(m) (m) (m) (2 LOPs) 
НТ. 332 от 11.8 1-1 * 

1.2 bus 3.0 х 
Т. 35220001 ж TT * 

I DN 1 
2 2:2 332 03 2.2 0.8 х 

| 1.2 me" 2.1 х 

ЯГ: онд: 4333 0.2 3.5 [NN : 

0.2 5.0 2.4 x 
ЖЕРІ 2 332 0-04 51 0.8 * 

още 4.9 1.9 x 
вида 333 0207 212 ШІП х 

ORI Ди 286 х 
[mel Ae 37 50 5.0 5.0 90° 
ime lee 353 50 5.0 5.0 909 


* NOTE: For positions which could be computed by more 
than one combination of LOPs (e.g., 4 combinations of 
3-range LOP fixes exist when 4 ranges have been 
measured), the upper line for each day shows the 
combination of LOPs giving the minimum mean gp's 

and the lower line for each day shows the —— 
of LOPs giving the maximum mean gp's. 


VE 


1.5 m on day 332-апа421:9 тог 21127140 20200200 
Observing both a range and an azimuth, which intersect at 
90°, at each of the two stations causes tEhe бр values tombe 
more consistent. When the two normal LOPs from one station 
are combined with two normal LOPs from another station a 
strong position is determined. Since an angle of 
intersection of 90° is the best possible шог sy, for a 
two-LOP fix, then this position is made from two accurate 
sets of LOBs. 
с.  Two-Range and One-Azimuth LOPs 

Overall mean values of Op of 1.4 to 2.0 m were 
determined when two-range and one-azimuth LOPs were used for 
positioning.  Variability exists in the range of gp Values 
observed due to the geometry problems stated previously. 
Two of the LOPs intersect at 90? and then a third LOP is 
added to this accurate configuration to determine the 
СПпиеес- БОР ме. 

4. One-Range and One-Azimuth LOPs 

The observed бр for the two-LOP Range and 
Azimuth method varied very little due to the minute change 
in the intersection angle of the LOPS. Since the ranges 
were measured from eccentric stations for three of the four 
control sites, the angle of intersection was slightly 
different than 90°. A value of 5 m for the standard error 


of position was determined for most of the position fixes. 
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Due to the consistency of the intersection angle 
of the LOPs, the range-azimuth LOPs method is prefered for 
large-scale hydrography. With the present interpretation of 
the IHO specifications, the 5-m бр would be acceptable for 


positional accuracy. 


C. COST COMPARISONS OF POSITIONING METHODS 

The increased cost of using MLOP may outweigh the 
ДЯ of an increase in accuracy. Conventional two-LOP 
methods require the minimum equipment and personnel to 
produce a hydrographic survey. The continuing increase in 
marine positioning requirements may force the hydrographer 
of the future to decide what price should be paid to produce 
a more accurate survey. 

Component costs of the equipment used in this thesis 
were obtained in June 1985 from manufacturer representatives 
(Table VIII). Prices quoted are the government contract 
rate in U.S. dollars. Table IX combines the costs of all 
the components depending on the method used to position the 
vessel. Only the four-LOP and two-LOP methods.were 
addressed as the costs of the other methods can be inferred. 
The costs for the two-LOP method Въ the data will be 
acquired by the conventional nonautomated procedure commonly 
used on large-scale surveys. Equipment which is common to 
all methods (i.e., vessels, tripods, signals, etc.) have 


been ignored in the tabulation. 
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TABLE VIII 


COMPONENT EQUIPMENT COSTS 


COMPONENT COST($) 
Del Norte Trisponder: 
Model 542 DDMU 14,500 
Master 9,500 
Remote 9770100 
RACAL-DECCA AUTOCARTA ІІ 24,750 
Houston Instrument Plotter 6,500 
Wild T-2 Theodolite 6,500 
Portable VHF Radio 800 


The price of оси т the four-range MLOP automated 
method approaches $100,000. In comparison, equipment for 
the four-azimuth MLOP method is about one-third the price of 
the four-range MLOP method. The difference in equipment 
costs in the two methods is very deceptive. To acquire four 
azimuths, four observers are needed to track the vessel. 

The cost for this additional personnel over the automated 
ranging system could approach the difference in equipment 


costs in just over one year. The labor costs involved would 
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TABLE IX 


COST COMPARISON OF EQUIPHENT 


Four-Range MLOP 


540 DDMU 
Master 
Remotes (4 
AUTOCGARTA «1 
Plotter 


@ $9,600) 


TOTAL 
Four-Azimuth MLOP 


T-2 Theodolites v 4: 500) 
Portable Radios (4 00) 


520-4200 


p= 


TOTAL 
Two-Range & Two-Azimuth MLOP 
540 DDMU 
Master 
Remotes -— с 99, 600) 
UCT 
-2 Theodolites 6 ë d: 500) 
Portable Radios 00) 
TOTAL 


Two-Range LOPs 


540 DDMU 
Master 
Remotes (2 
АОТОСАКТА 1 


С $9,600) 
TOTAL 
Two-Azimuth LOPs 


T-2 Theodolites 6 500) 
Portable Radios 800 0) 
TOTAL 


Range and Azimuth LOPs 


540 DDMU 
Master 

Remote | 
T-2 Theodolite 
Portable Radio 


TOTAL 
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depend upon the skill and experience of the personnel used 
to acquire the data. For four observers paid $8.00 per 
hour, for an 8-hour day, for 260 days per year, the increase 
due to labor would be 566,560 (4 x 8 x $8 x 260) per year. 
Added to the equipment costs of $29,200, the total price 
exceeds $95,000. 

Another expense associated with the four-azimuth MLOP 
method would be the cost of the hardware needed to process 
the data using a least squares algorithm. A data processing 
computer would be necessary to perform the repeated 
iterations required by the algorithm. Additionally, a major 
disadvantage of the nonautomated data acquisition method is 
that observed azimuths must be hand logged into the data 
processing system. This manual intervention increases the 
existance of numerous blunders and involves additional labor 
costs which could approach the cost of acquiring the 
Еопг-аг1 ша LOPS. 

With these facts in mind, the overall costs of the 
four-azimuth method could approach the cost of the automated 
four-range method. Given this circumstance, the optimal 
system would be to use an automated ranging method to 
increase production, reduce blunders in recording the 
observations, and have better control of the survey 


operations. 
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DR OPTIMIZED SYSTEM 

One of the objectives of this thesis was to determine an 
optimal system for the acquisition of positional data in 
hydrography. Results from the data analysis proved 
conclusively that MLOP do improved fix accuracy over two 
LOPs. Every effort was made to duplicate the actual field 
work currently done by NOS so that a valid comparison of 
methods could be made. 

To adequately discuss an optimal system a few 
considerations must be made. First, a consideration of the 
accuracy requirements must be made. This accuracy 
БИЕ semen: will be dependent upon the purpose of the 
Survey. Even with the initial accuracy estimates used for 
this thesis, all of the MLOP and the range-azimuth two-LOP 
methods will meet the 5-m positional accuracy. Therefore, 
the optimal system would be dependent upon the final product 
the hydrographer wishes to produce. If automation is the 
overriding criterion, then a ranging MLOP method wouid 
ће fulfill the needs of the survey. An automated 
system, — as AUTOCARTA II, can be bought or leased 
off-the-shelf, to provide the user with an automated MLOP 
system. Automation frees personnel to do other tasks since 
the need for shore observers is eliminated (except for 
station maintenance). By acquiring the data automatically 


on board the survey vessel, the existance of blunders are 
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minimized, thereby decreasing the processing time. Since 
the capability exists for pre-selecting track lines and 
controlling the movements of the vessel, more production is 
possible with an automated system. Overall the use of range 
MLOP and an automated DAS is the most feasible method to 
conduct accurate near-shore hydrography. 

Second, if the purpose of the survey dictated that the 
accuracy of position must be to within 1 m, or less, then 
more effort must be made to acquire this accuracy. Ewell 
when a 4-m standard error was used for each LOP, this thesis 
found that the Op for 63% of the positions determined by 
four- and three-azimuth LOPs were less than or equal to 1.0 
m. Although the method is susceptible to numerous errors, 
is labor intensive, and logistically complex, one could 
produce a survey of a mearshore area to within 1.0-m 
accuracy with good field procedures. 

Third, cost must be considered. The benefit of using 
Just two LOPs is that minimal equipment and labor costs are 
involved to acquire adequate data. The use of MLOP will 
inczdiili the initial cost of equipment at least twofold. 
The cost may be prorated over time so in 4 to 5 years the 
benefit of producing a higher quality product would 
compensate the hardware costs. Additionally, the cost to 
establish more shore control must be considered. More LOPs 


would mean more control stations are required than presently 
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needed for conventional hydrographic surveys. Detailed 
survey planning could minimize the cost by choosing the 
optimum number of stations needed for the method of data 


acquisition. 
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VI. CONCLUSIONS 


A. MLOP FOR ON-LINE CALIBRATION OF RANGING SYSTEMS 
Short-range microwave survey systems, such as the Del 
Norte Trisponder, require a static base-line calibration 
over a known geodetic distance to eliminate systematic 
errors. Daily system checks must also be made during the 
course of the survey to detect any calibration drifts over 
time. Daily system checks can be made in either a static or 
dynamic mode. Static checks can be made by taking fixed 
point observations alongside a known point, three-point 
sextant fixes with a check angle, theodolite intersections 
from shore control, or range-azimuth positions using a total 
Station EDM unit. Dynamic checks can include steering a 
range with predetermined crossing angles or underway 
Simultaneous theodolite intersections [Holder, 1983]. 
Experience has shown that whatever system is used to 
perform daily calibration checks of the equipment, a 
substantial amount of time is consumed on this exercise in 
quality control. Weather permitting, daily system checks 
are made prior to and after each day's hydrography. Under 
ideal conditions, the amount of time consumed to obtain the 


required calibrations usually is l1 hour of the survey day. 
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равен пие таа сарено Е tusine MLOP for position control 
in hydrography is that redundant data are acquired 
Simultaneously on each fix. By having a third or fourth 
LOP, errors in position could be detected more readily than 
in the conventional survey mode where everything was assumed 
to be correct until it was too late. Thus the need to halt 
Survey operations to make a daily system check could be 
eliminated. Suppose four-range MLOP were to be used to 
position the vessel during a survey. Once the system was 
iine ana warmed-up, a number of positions could be 
obtained by the automatic system while underway to the 
working grounds. If a large standard error of position, 
Op» was determined by the software within the computer, 
then an examination of the аана = for each LOP could be 
examined to determine which range or ranges were in error. 
This error could be from equipment malfunctions, errors in 
the locating of shore stations, or erroneous ranges due to 
low signal strength, reflections, etc. Additional fixes 
could ASR um. or reject the existance of a substantial 
error. Once isolated, the unit (or units) causing the error 
can be removed or replaced. Here the advantages of having 
four-range MLOP as control can be realized. Should one of 
the units fail, there would still be a redundant range 
acquired allowing for a least squares computation. If for 
some reason two units failed, operations could continue 


uSing conventional two-range LOPs. 
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B. FUTURE ASPECTS AND CONSIDERATIONS 

The use of MLOP for hydrographic sounding line control 
should be the primary method for future nearshore surveys. 
As advances in electronic technology continue, the economy 
of miniature computer installations will make automated 
hydrographic data acquisition systems more feasible. NOS is 
presently in the process of revitalizing the outdated 
HYDROPLOT/HYDRLOG systems onboard its hydrographic vessels. 
The Shipboard Data — III (SDS III) is currently bes йн 
built to increase automation and acquisition of hydrographic 
data. One of the capabilities of SDS III will be to acquire 
and process MLOP for sounding line positions. 

The methods presented in this thesis for data 
acquisition are based on the use of components, such as Del 
Norte Trisponders, which have been in existence since the 
early 1970's. This thesis has demonstrated that using these 
systems in a MLOP configuration will substantially increase 
the accuracy of the survey. However, other systems exist 
which should be considered in conjunction with conventional 
УДОВА 

Krupp-Atlas Elektronik introduced a new short-range 
dynamic positioning system, Polarfix, in 1982. The 
advertised [Wentzell, 1983] advantages of the system are: 

(a) Range-Azimuth fixing to 0.1l-m accuracy. 
(b) Automatic operator-free tracking of the vessel. 


(c) Portable shore stations. 
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Since the Polarfix system has the potential to alleviate 
some of the accuracy problems associated with the present 
two-LOP methods, additional studies should be made to 
ascertain the manufacturer's claim as to the proported 
accuracy of the system. 

Differential GPS is a method which should be considered 
for large-scale hydrographic surveys. The theory has proven 
GPS extremely accurate (within centimeters) for control 
Sites ashore where receivers acquire positional data from 
numerous passes of the satellites.  Accuracies on the order 
of 2 to 3m may be attained underway once the theory and 
technology have been developed. | 

Future positioning systems may be integrated so that a 
combination of methods and systems can be analyzed to 


determine the "best" fix. 
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APPENDIX A 


FORTRAN PROGRAMS FOR DATA PROCESSING 


The following listed programs were used to process the 
LOP data for this thesis. Output files were generated to 
list the standard error of the adjusted coordinates, Te and 


бу; the correlation coefficient, the semi-major and 


Oxy? 
зеш?-тш1 пот” ажез”о# "е егготте! зе о eM NINE 
orientation of the error ellipse, §; and the standard error 


of position, Ср: 


А22Т.ОР 
Program AZ2LOP computes the positional accuracy of a fix 


obtained from two azimuth LOPs. This program is a 


modificaton of program UCOMPS. 


-LSQAZ4 

Program LSQAZ4 computes positional accuracy using four 
azimuth LOPs. This program is a modification of program 
SILVAl [Sliva, 1982]. The modification permits handling 


large amounts of data. 


RAZ2LOP 

Program RAZ2LOP computes positional accuracy of a fix 
obtained from one range LOP and one azimuth LOP. This 
program is a modificaton of program UCOMPS. It includes the 
error propagation algorithm. 
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RR+AZ 

Program RR*AZ computes the positional accuracy of a fix 
obtained from two range LOPs and one azimuth LOP. This 
program is a modification of program SILVA3 [Sliva, 1982]. 


The modification permits handling large amounts of data. 


RR2LOP 
Program RR2LOP computes positional accuracy using two 
range LOPs. This program is a modification of program 


UCOMPS. 


R2+AZ2 

Program R2+AZ2 computes the positional accuracy of a f1x 
obtained from two range LOPs and two azimuth LOPs. This 
program is a modification of Soon ои f Silva) 1982]. 


The modification permits acceptance of another azimuth LOP. 


R3+AZ 

Program R3+AZ computes the positional accuracy of a fix 
obtained from three range LOPs and one azimuth LOP. This 
расова о а modification of program SILVA3 [Silva, 1982] 
which computed a least Squares position from two range LOPs 


and one azimuth LOP. 


SLEAST 
Program SLEAST (Appendix B) computes the least squares 


position and positional accuracy of a fix from three or four 
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range LOPs. The program was written by Lt. Nick Perugini 


and the author from the algorithm presented by Cross [1981]. 


UCOMPS 

Program UCOMPS computes coordinates of a point in either 
Geographical Position (G.P.) or X and Y grid coordinates 
from range-range, range-azimuth, or azimuth-azimuth LOPs. 
It is a utility package program used by Hydrographic 


Sciences students at NPS. 
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APPENDIX C 


FOUR-RANGE MLOP FIX ACCURACIES 


EXAMPLE DATA OUTPUT FROM PROGRAM SLEAST 
POSITION FROM FOUR RANGES 


CONTROL STATIONS WERE SQUARE, CONK, USE MON, GEOCEIVER 


FIX SIG X SIG Y SIG XY SIG A SIG B THETA SIGMA P 


1 0.90 1.3 0458 eee 0.81 20200876 
2 аза 10 ә 101, -20 ТОШ 
3 1.39 1.62 000.76 ПОЗНЕ 55-250 МИ 
4 1.56 1.79 0785 198 Ч? -315 088 
9 0.32 0.39 90 029 901 200 tome 
10. 0.76 0.97 Оо (87 70. 0.30 1.23 
11 40559 0.76 -0027 о 18: 2:36 0.9 
137 0.77 0.95 40.04 о 5.69 1.19 
14 0.59 0.80 -0.04 0.80 0.59 7.42 0.99 
15 10.597 70.80 2-20 7000 О: 9.04 1.00 
17 0.32 0.45 -0.02 Оо ОТ 10949421 
18 0.47 0.67 -0.06 069 0296 13.39 ОЙ” 
20 0.4 1.39 -ол 15790959 5 15 
21 1.09 1.25 -0.69 1.67 95 30 ҰС 
22 1.15 1.28 -0.64 1.63 10Іі 5502: ТБ 
23 1.70 1.84 -1.14 2.23 151 5871 204 
24 1.69 1.79 -0.92 2.09 1.52 39.40 2.46 
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APPENDIX D 


FOUR-AZIMUTH MLOP FIX ACCURACIES 


DATA PROCESSED USING A MODIFICATION OF PROGRAM SILVAL 
FROM SILVA (1982) TO HANDLE LARGE DATA FILES. 

EXAMPLE DATA OUTPUT FROM PROGRAM LSQAZA. 

POSITION FROM FOUR AZIMUTHS. 


CONTROL STATIONS WERE SQUARE, CONK, USE MON, GEOCEIVER. 


FIX# SIG X SIG Y SIG XY SIG A SIG B THETA SIGMA P 


fo 008 0.15 -0.00 0.16 0.08 97.59 0.17 
77° 0.04 0.07 0.00 0.07 0.04 . 90.12 0.08 
тэг 0597 005 992 001 83.75 1.70 
79 0.31 0.41 0.02 0.41 0.30 Тоја ОТУ 
80: 0.23 0.28 0.01 0.29 0.23 1582798 097 
81 0.80 0.88 0.05 0.89 0.79 71.48 1.19 
82 0.43 0.45 0.00 0.45 0.43 76.41 0.62 
ЕП ЖЕП 73 001 0:3 0.42 129.29 0.60 
оо -002 053 0:9 13763 0.72 
О О 001 5027 027 15585 0.36 
ЕРЕ 259033 002 1355 0.30 132.84 0.47 
87 0.39 0.40 -0.02 0.42 0.37 131.27 0.56 
88 0.56 0.55 -0.03 0.59 0.53 135.74 0.78 
89 0.43 0.40 -0.01 0.44 0.39 158.24 0.59 
ЕЖ 902-009 122 90.99 170.41 1.57 
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APPENDIX F 


TWO-RANGE AND TWO-AZIMUTH MLOP FIX ACCURACTES 
EXAMPLE DATA OUTPUT FROM PROGRAM R2+AZ2 


POSITION FROM 2 RANGES AND 2 AZIMUTHS 


CONTROL STATIONS WERE SQUARE, SQUARE ECC., USE MON, USE MON 


FIX# SIG X SIG Y SIG XY SIG A SIG B THETA SIGMA P 


1 1.47: 166 Ой 2 71272. 68.63 2.22 
2 1.81 2.04 072 47054 ММ 67.92 2.73. 
G 1.97 2.16 0.45 2.22 1.93 65.66 2.93 
6 0.97 1.4 omo "027 85 63.17 1.43 
7 0.84 0.89 005 9907 Шо 61.98 127 
8 0.57 0.60 902 001 111 60.74 0.82 
9 0.85 0.88 0.05. 0.90 0.84 59.44 1.23 
10 0.99 1.01 0.05 1.03 0.97 58.19 1221 
11 0.71 072 002 0 70 56.92 10 
13 0.71 0.71 0.00 5071 Ол 147351 198 
14 0.49 0.40 -0.00 0.41 09627 2002 700 си 
15 0.83 0.32 -0.03 0.85 М: Cico ШЫ 
17 0.58 0.57 -0.03 0.60 0.55 138.59 078 
18 0.55 0.55 -0.04 0.58 0.52 137.41 о 
20 0.24 0.24 -0.01 0.26 0.22 13506 OME 
21 1.00 1.11 -018 1.17 0999 527 
22 На 0.84 118.64 1.30 
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